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Abstract. In this article, we give sharp bounds for the Euler- and trapezoidal dis- 
cretization of the Levy area associated to a d-dimensional fractional Brownian motion. 
We show that there are three different regimes for the exact root mean-square conver- 
■ gence rate of the Euler scheme. For H < 3/4 the exact convergence rate is 

CS| ! where n denotes the number of the discretization subintervals, while for H = 3/4 it is 

• n~^{log{n)y^'^ and for H > 3/4 the exact rate is n~^. Moreover, the trapezoidal scheme 

D , has exact convergence rate 7j-2^f-(-i/2 fgj. jj > 1^2. Finally, we also derive the asymptotic 

error distribution of the Euler scheme. For H < 3/4 one obtains a Gaussian limit, while 
^T) , for > 3/4 the limit distribution is of Rosenblatt type. 



1. Introduction and Main Results 

Let B = {B^^\ . . . , B^^^) be a rf-dimensional fractional Brownian motion (fBm) with 
Hurst parameter H G (1/4, 1) indexed by M, i.e. B is composed of d independent centered 
Gaussian processes whose covariance function is given by 

> : Rh{s, t) = l + -\t- .r^) , s,teR. 

: 

0> . For an arbitrary T > 0, a typical differential equation on [0, T] driven by B can be written 

O ■ rt 

c4- Yt = a+ a{Y,)dB,, t e [0,T], (1) 

o ; Jo 

^ ■ where a G M*^ is a given initial condition and a : M" ^ M"''^ is sufficiently smooth. During 

the last years, the rough paths theory has allowed to handle several aspects of differential 
equations like ([T]), ranging from existence and uniqueness results (see [7l[l4] for equations 
^ I of type dl]) and [HI [IHl [19] for extensions to other kind of systems) to density estimates [4] 

■ or ergodic theorems [TT] . 

It is also important, and in fact at the very core of the rough path analysis, to derive 
good numerical approximations for fractional differential equations like ^ . This problem 
has so far been considered in three type of situations: (i) When > 1/2, it is proved 
independently in [6] and [H] that the Euler scheme associated to equation ([Tl), based 
on a grid {iT/n; i < n}, converges with the rate for arbitrarily small e > 0. 

The exact rate of convergence of the Euler scheme is computed in [18] in the particular 
case of a one-dimensional equation, (ii) In the Brownian case H = 1/2, there exists 
a huge amount of literature on approximation schemes for SDEs, and we just send the 
interested reader to the references [13 [T5| for a complete overview of the topic, (in) For 
1/3 < H < 1/2, the rough path strategy in order to solve equation (P), see e.g. [71 l8l fT4]. 
tells us that one should use at least a Milstein-type scheme in order to approximate its 
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solution. Moreover, it can be easily seen that for < 1/2 the standard Euler scheme does 
not converge and in fact explodes for stepsizes going to zero, even in the one-dimensional 
case. Indeed, consider for instance the one-dimensional SDE 

dXt = Xt dBt, Xq = 1, 

whose exact solution is given by Xt = exp{Bt). The Euler approximation for this equation 
at t = 1 is given by 

n-l 
A:=0 

So for n G N sufficiently large and using a Taylor expansion, we have 

n-l 

X, - Xf"^ = exp(Si) - exp ( 5^ log(l + {B^k+i)/n - %n))) 

fc=0 

n-l 



exp(Si) - exp (^1 - ^ 5Z \^{k+l)/n - Bk/n\^ + Pn) : 

k=0 



Prob 

where pn — > for n ^ oo for > 1/3. Now it is well known that 



|-D(fc+l)/n — -Dfc/nl ^ OO 



fc=0 



ior H < 1/2, so we have x[^^ oo. However, Milstein-type schemes are known to be 
convergent for such a one-dimensional equation, see [9]. 

For the general multi-dimensional equations of type ([T]), a Milstein-type scheme is 
studied in set Yq = a, and for a grid given by tk = kT/n, k = 0, . . . ,n — 1, let 

n 

Yt,,, =yt, + E - Bt,) (2) 

i=l 

i,j=i •'^k 

for /c = 0, . . . , n — 1, where V^^^ is the differential operator Yl'i=i ^f^ ^^^i ■ Davie then proves 
that this scheme has convergence rate fi-i'^H-'^)+<^ ^ and this result has been extended in [7] 
in an abstract setting, to higher order schemes for a rough path with a given regularity. 

The above Milstein-type scheme ([2]) requires knowledge of the iterated integrals 

XP= fi^Wrffip), tG[0,T], z,j = l,...n, (3) 

whose explicit distribution is unknown for i ^ j. Thus discretization procedures for 
([3]) are crucial for an implementation of this numerical method. This has already been 
addressed in [5], where dyadic linear approximations of the fBm B are used in order to 
define a Wong-Zakai-type approximation X"^ of X. In the last reference, the process X^ 
is shown to converge almost surely in p-variation distance, and the (non-optimal) error 
bound 

E|X^ - Xt\^ < C ■ 2-"(4J^-i)/2 
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is also determined. The current article takes up this kind of program, and we consider 
the approximation of 



^0 



(4) 



by the Euler- and a trapezoidal scheme based on equidistant discretizations. 

For the approximation of ([4]) the standard Euler method has the explicit expression 

n-l 
i=0 

The results we obtain for the Euler scheme are then of two kinds. First we determine the 
exact L^-convergence rate. 

Theorem 1.1. Let Xt defined by ^ and its Euler approximation XJj^ given by expression 
Define the constants aj{H), for j = 1,2,3 by 

a,{H) = - ( (3i2H, 2H) + + - ( (1 - 2^^) + + ^^"l 

^ 2V ^ AH-lJ 2y ' AH-l AH-l) 

+ H [\y'^\l+yr^'-y'^^'\l + yndy 
Jo 



and 



am = am + ^!M^C(4 - 4i7), a.^H) - ' ^'^'^ 



2 " ^ 4 4:H-3 
Then we have 

ai(/7)-T4^-n-^^+i + o(n"4^+i) for H e (1/4,1/2), 

_ 2_ J a2(i^)-T4^-n-4^+^ + o(n-4^+i) for H e (1/2,3/4), 

£.|At a^i - <; _|_ . . iQg(^)^-2 ^ o(log(n)r2-2) for H = 3/4, 

as{H) ■ T"^" ■ n-^ + ofn'^) for H e (3/4,1). 

Observe that for the case H = 1/2, i.e. for the approximation of the Wiener Levy 
area, one obtains by straightforward computations that E|Xt — = ^ ■ n~^, which 
is compatible with our Theorem 11.11 since 

lim ai(H) = lim a2(H) = -. 

H^l/2,H<l/2 H^l/2,H>l/2 2 

The convergence rate breaks up into several regimes which are reminiscent of the cases 
obtained in [2H[24] concerning weighted quadratic variations of the one-dimensional fBm. 
In particular, the convergence rate does not improve for H > 3/4, i.e. is equal to 
independently of H. Finally, note that our study starts obviously at H = 1/4^, since the 
Levy area is not even defined for if < 1/4. 

Using a trapezoidal rule for the approximation of the integral leads to the following 
scheme, which coincides with the Wong-Zakai approximation used in [5]: 

^ n— 1 

^T= {^iT/n + B{ill)T/n) (^(?il)T/n " ^iT/n) ' (6) 

i=0 

This trapezoidal scheme avoids the "breakdown" of the convergence rate of the Euler 
scheme for if > 3/4. 
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Theorem 1.2. Let H > 1/2. Then we have 
where 

Note that the constant a^i^H) could also be expressed in terms similar to ai{H). How- 
ever, we think that this gives no further insight and thus we omit it here. We strongly 
suspect that the root mean square convergence rate n"^^"*"^/^, which is obtained by this 
trapezoidal scheme, is the best possible. In other words, we conjecture that the conditional 
expectation of Xt given -Br/m B2T/n, ■ ■ ■ Bt satisfies 

E\Xt - E{Xt I Br/n, B2T/n, ...Bt)\' = C{H) ■ T^^ ■ n'^^+i + o(n-^^+i), 
where C{H) > 0. 

The third result in this article is a refinement of Theorem 11.11 meaning that we obtain 
a limit theorem for the asymptotic error distribution of the Euler scheme. 

Theorem 1.3. Let X^, XJ}, and ai{H) , a2{H) , a^{H) defined as above. Moreover, let Z 
he a standard normal random variable. Then : 

(1) Case 1/4 < H < 3/4.- the following central limit theorems hold: 

1- „2H^i/2(^ ^n. (D r ^/MH)T'''-Z for He (1/4,1/2), 
™^ (^T-X^) - I He {1/2,3/4) 



and 



hm n(logH)-V2(x^_x?) -^T'^ ■ Z 



for if = 3/4. 

(2) Case H > 3/4.- let Ri and R2 be two independent Rosenblatt processes (see Sec- 
tion 5 for a definition). Then it holds 

lim n {Xt - X^) ^2ai{H)T^" ■ {R^ - R^). 

n— »oo 

Let us say a few words about the methodology we have adopted in order to prove Theo- 
rem [HSl It should be mentioned first that we have used the analytic approximations intro- 
duced in [25] in order to define the Levy area X, which allows to use some elegant complex 
analysis methods for moments estimates in this context. Then, for H e (1/4,3/4), the 
central limit type results are obtained through the criterion introduced in [22] for random 
variables in a fixed chaos. For this we control the fourth moments of X with the help 
of (Feynman) diagrams. For the case H > 3/4 we proceed in a different way. Here the 
Milstein approximation of X^ performs better than the Euler method. Then expressing 
the differences between both schemes as the sum of quadratic variations for two inde- 
pendent one-dimensional fBms, thanks to a simple geometrical trick given in [20], one 
obtains the limit theorems for if > 3/4 using the results of [24]. In particular, this leads 
to the Rosenblatt type limit distribution as in [24]. For the trapezoidal scheme, whose 
error behaves like the second order quadratic variations of fBm, see e.g. [2], a central 
limit theorem could be also derived using the criterion in [22], but we omit this here for 
the sake of conciseness. 
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The remainder of this article is structured as follows. Integrals with respect to the 
fractional Brownian motion will always be understood as limits of analytic integrals as in 
[25] . We thus recall the definition of the analytic fBm, as well as some preliminaries at 
Section [21 Section [3] contains the proofs of Theorem 11.11 and II. 2[ The proof of Theorem 
11.31 is given in Sections [4] and El 

2. Definition of the analytic fbm and preliminaries 

This section is devoted to recall the definition of the fractional Brownian motion in- 
troduced in [25], and to state some of the properties of this process which will be used 
in the sequel. All the random variables introduced here will be defined on a complete 
probability space P), without any further mention (notice the unusual notation 

lA for our probability space, due to the fact that the letter VL will serve for the complex 
domains we consider in the sequel). The following kernels will also be essential for our 
future computations: 

Definition 2.1 (?7-regularized power functions). For /? G M \ Z and r] > let 

[x]^'' = {±ix + rif and = 25R[x]±'^ = [a;]+'^ + [a;]-'^ 
Then, for 7] > and x,y eR, define K''^{r]] x, y) as 

K'Hv, X, y) = ^J^"'^^ (±i(x - ,) + ,)^^- = ^^^[x - y]^"-'. 

2 cos IT H 2 cos IT H ' 

Set also 

K'{r];x,y) := 2^K''^{r];x,y) = K''^{r];x,y) + K'~{r];x,y). 
Notice that the above kernels are well-defined on our prescribed domain M!J_ x M x M. 

2.1. Definition of the analytic fBm. The article [25] introduces the fractional Brow- 
nian motion as the real part of the trace on M of an analytic process T (called: analytic 
fractional Brownian motion [23]) defined on the complex upper-half plane 11+ = {z E 
C; '^{z) > 0}. This is achieved by first noticing that the kernel K'{ri) is positive definite 
and represents (for every fixed 77 > 0) the covariance of of a real-analytic centered Gauss- 
ian process with real time-parameter t. The easiest way to see it is to make use of the 
following explicit series expansion: for > and z G n"*", set 

n-, [WZM l n2-2H + k) fz + iV^-Wz-iy 
f'^'^-^ V 2cosvri/ V T{2-2H)k\ J ' 

where F stands for the usual Gamma function. Then these functions are well-defined on 
n+, and it can be checked that one has 

^fk + fk + = (^(^i + ^2);a;,?/j . 

fc>0 ^ ^ 

Define more generally a Gaussian process with time parameter z G n"*" as follows: 

T\z) = Y,h{z)^k (8) 

fc>0 

where (^fc)fc>o are independent standard complex Gaussian variables, i.e. E[^j^fc] = 0, 
E[^j^fc] = 5j^k- The Cayley transform z maps n"*" to where V stands for the 

unit disk of the complex plane. This allows to prove trivially that the series defining V 
is a random entire series which may be shown to be analytic on the unit disk. Hence the 
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process T' is analytic on 11+ . Furthermore note that, restricting to the horizontal line 
M + i^, the following identity holds true: 

E[r'(x + ir//2)r% + V2)] = K'^-{ri; x, y). 

One may now integrate the process F' over any path 7 : (0, 1) — > 11+ with endpoints 
7(0) = and 7(1) = z E 11+ U M (the result does not depend on the particular path 
but only on the endpoint z). The result is a process F which is still analytic on 11+ . 
Furthermore, one may retrieve the fractional Brownian motion by considering the real part 
of the boundary value of F on M. Another way to look at it is to define Ttit]) := F(t + i?]) 
as a regular process living on M, and to remark that the real part of F(?7) converges when 
77 ^ to fBm. In the following Proposition, we give precise statements which summarize 
what has been said up to now: 

Proposition 2.2 (see [25l[23]). Let T' be the process defined on 11+ by relation 

(1) Let 7 : (0, 1) — i> n+ be a continuous path with endpoints 7(0) = and 7(1) = z, 
and set = J^T'^du. Then F is an analytic process on 11+. Furthermore, as 
z runs along any path in 11+ going to t G M, the random variables F^ converge 
almost surely to a random variable called again Tf. 

(2) The family {Tf, t G M} defines a centered Gaussian complex-valued process whose 
paths are almost surely K-Holder continuous for any k < H. Its real part Bt : = 
29fJFt has the same law as fBm. 

(3) The family of centered Gaussian real-valued processes Btirj) := 23?Ft+i^ converges 
a.s. to Bt in a-Holder norm for any a < H, on any interval of the form [0, T] for 
an arbitrary constant T > 0. Its infinitesimal covariance kernel 'EiB'^{7])B'y{rj) is 
K'{r];x,y). 

2.2. Definition of the Levy area. Let us describe a natural possible definition of the 
Levy area associated to F. Since the process Bt{r]) := 23fJFt+i^ is a smooth one, one can 
define the following integral in the Riemann sense for all < s < t and 77 > 0: 




It turns out that A{r]) converges in some Holder spaces, in a sense which can be specified 
as follows. Let T be an arbitrary positive constant, Cj be the set of continuous complex- 
valued functions defined on [0,T]-', and for /i > 0, define a space of /x-Holder functions 
on [0,r]2 by 

\\f\\,:= sup and C^{V) = {f e C2{n;Vy, \\f\\, < c^} . (10) 

The /^-Holder semi-norm for a function g E Ci is then defined by setting hst = gt — 9s as 
an element of C2, and \\g\\^ := \\h\\^ in the sense given by ffTOll . 

According to [25], [23j, the Levy area A oi B can then be defined in the following way: 

Proposition 2.3. Let T > be an arbitrary constant, and for s,t E [0, T]"^, r] > 0, define 
Astiv) C"^ equation Consider also < 7 < if. Then: 

(1) For any p> 1, the couple (5(r/), ^(r/)) converges in Lp{Q;CJ{[0,T];R)xC2^{[0,T]'^;R)) 
to a couple {B,A), where B is a fractional Brownian motion. 
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(2) The increment A satisfies the following algebraic relation: 

At - A„ - At = ^ - (fii^^ - Bi'^) , 

for s,u,t E [0, T] . 

Notice that the algebraic property (2) in Proposition 12.31 is the one which qualifies A 
to be a reasonable definition of the Levy area of B. 

It will be essential for us to estimate the moments of A. For this we will use the 
following definition: 

Definition 2.4. For rj > and 01,02 G M, let us define the function Kai^a2{V'r: ') on 
RxR by 

Ka^^a2iV:XuX2) = dyi dy2K'{ri;yi,y2). (11) 

J ai J Ci2 

Notice then that, invoking the conventions of Definition \2.1[ we have 

+ [0i-02]f). (12) 

We also state the classical Wick lemma for further use. 
Proposition 2.5. Let Z = {Zi, . . . , Z2n) be a centered Gaussian vector. Then 

N 

E[Z^--.Z2n]= E 1[E[Z,,^Z,,^^,] (13) 

(n,fi2),...,(i2]v-i,«2jv) i=i 

where the sum ranges over the {2N — 1)!! = 1-3-5-- - {2N — 1) couplings of the indices 
1,...,2N. 

We can now give the announced expression for the moments of A{ri) (recall that A{ri) 
is defined by Q): 

Lemma 2.6. Let N > 1 and {sj, tf, i < 2N} be a family of real numbers satisfying Si <ti. 
Then 

dxi - ■■ / dx2N (14) 

-1 ■^S2JV 

N N 

(n,i2),.--,(i2iV-l,«2Jv) O'l J2),.--,(j2]V-l,j2Jv) fc=l k=l 



E 



2N 



Proof. By definition of the approximation A{r]), we have 



E 



'2N 



2N 

n 



dyjE 



B'}P{V)B'^^^\V) ■ ■ ■i?:S(^)^;S(^)l (15) 



Our claim stems then from a direct application of Proposition l2.2l point (3), Proposition l2.5l 
and Definition 12. 4[ 

□ 
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2.3. Analytic preliminaries. We gather here some elementary integral estimates which 
turn out to be essential for our computations. The first one concerns the behavior of the 
kernel -ft'ai,a2 given at Definition 12. II when |ai — |a2 — X2I are of order 1 and |xi — X2I 
is large. 

Lemma 2.7. Assume 77, |xi — ai|, \x2 — 02! < 1 and \xi — X2\, |ai — X2I, |a2 — \ai — 02] 
are bounded from below by a positive constant C. Then 

|/Cj,„2(?7;xi,X2)| < C (min(|xi - X2I, |ai - X2I, |a2 - ^i], |ai - 02!))^'^"^ . 

Proof. The proof is elementary using the integral expression ( fT2l) for -ft'ai,a2- 

□ 

We shall also need to estimate convolution integrals of the form K'{rj] z, u)f{u) du or 
/o Kafi{ri; u)f{u) du. The following lemma gives a precise answer when / is analytic on 
a neighborhood of (a, h) for two given constants a, 6 G M, and multivalued with a power 
behavior near a and h. 

Lemma 2.8. (see |26| i Fix two real constants a, h with a < b, and let f be a function in 
L^{[a, b], C). Define another function (f) by (f) : z t-^ Jl^{—i{z — u)Y{u — a)'^ f{u) du with 
7 > -1 and /3 + 7 e M \ Z. Then: 

(1) Assume f is analytic in a (complex) neighborhood of s G {a,b). Then has an 
analytic extension to a complex neighborhood of s. 

(2) Assume f is analytic in a complex neighborhood of a. Then may be written on 
a small enough neighborhood of a as the multivalued function 

^{z) = {z-aY+^^^F{z) + G{z) (16) 

where both F and G are analytic. 

(3) More precisely, the following continuity property holds: let Q be a complex neigh- 
borhood of [a,b] and e G (0, 1/2). // / is analytic on a relatively compact domain 
Vt containing the closure Vt of Vt, then extends analytically to the cut domain 
^cnt := fl\{{a+R-)U{b+R+)) and writes {z-a)'^+'^+^F{z)+G{z) on B{a,e{b-a)) 
(F,G analytic) with 

sup 101, sup sup |G'|<Csup|/| (17) 

n^ut\iB{0,eib-a))UB{b,e[b-a))) B{a,e(h-a)) B{a,e(b-a)) Q 

for some constant C which does not depend on f . 

Proof. Points (1) and (2) follow directly from [26j, Lemmas 3.2 and 3.3. Point (3) may 
be shown very easily by following the proof of the above two lemmas step by step and 
using the analyticity of /. Note that (under the hypotheses of (3)) is analytic on the 
larger domain ^2 \ ((a + M_) U (6 + M+)), but the method of contour deformation used in 
the proof gives a bound for 0(2) which goes to infinity when z comes closer and closer to 
the boundary of Vt (hence the need for the relatively compact inclusion of Vl into Vt). 

□ 

We shall also need the following elementary lemma. Here and later on, we will write 
X < ?/ for X, y G M, if there exists a constant C > such that x < C ■ y. 

Lemma 2.9. Let > — 1 and < a < 6 < 1. Then: 

[ |t-ar|t-6|^ rft < 1 + |a-6|"+^+^ (18) 
Jo 
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Proof. Let aa{b) = max(0,2a — b) and ab{a) = min(l,2fe — a). Split the above integral 
into Jq"^^^ + (^-) + + Jl^''^"^^ + Ja^[a)- show that the integral over each subinterval is 
< 1 + \a — b\°'^'^~^^ (by symmetry, it is sufficient to check this for the three first subintervals 
only). Now, a simple study of the function 1 1-^ \t — a\/\t — b\ shows that c < < C on 

[0,(Ta(6)], so 

/ \t - al^'lt - b\^ dt < (t - 6)°+^ rft < (6 - a)"+'3+i + (19) 

Jo Jo 

Ua + p + 1 <0, resp. a + + 1 > 0, then this is < (6-a)°+^+\ resp. < 1. On [aa{b),a], 
one has c < 4^ < C this time, so 

o—a ' 

f |t-a|°|t-6|'^ rft < (6-a)"+^+\ (20) 

Finally, jj^ \t — a|"|t — b\^ dt = ^^r(Q!+'a+2)^'' ~ a)""^^"'"^, where F is the Gamma fimction. 

□ 

3. Mean square error computations 

This section is devoted to prove Theorem 11.11 and Theorem 11.21 We will start with 
the error of the Euler scheme, the Milstein-type scheme will be considered later on. Note 
that we can decompose the error of the Euler scheme as Xt — = XliLi "Ji^^ where the 
random variables J" are defined by 

l'{i+l)T/n 

J- = / (Bi'^ - ) = A^,T)/n,i^+l)T/n, ^ = 0, . . . , n - 1, 

JiT/n 

where Ast is obtained as the L^-limit of Astiv) according to Proposition [231 In particular, 
E[|Xt — X^l'^] = X]j j JJ^], which means that we are first reduced to study the 
quantities 'E[J^J'!^] in terms of i,j and n. Towards this aim, one can first remark that, 
since fBm is self-similar and has stationary increments, we have 

E[Jf j;] = T^^n-^^E[i,J,], (21) 

with 

h = [ {Bi'^ - 5f )) = A^,^+u 2 = 0, . . . , n - 1, 
We now show how to handle those terms. 

3.1. Some moment estimates. The preliminary results we need in order to prove The- 
orem 11.11 are summarized in the following lemma: 

Lemma 3.1. Let Aoi = bP dB^^^ and Au = J^{bP - B[^^)dBP be the double 
iterated integrals with respect to B obtained by applying Proposition \2.3[ Define 

ci = E [lAiH , and C2 = E [AoiAu] ■ 

Then we have 

c, = ^((3{2H,2H) + ^^Y (22) 
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and 



1 



2H - 1 



H2 



AH 



H 

y 



Jo 



2H-1 _ y2-f/-l| 



y + indy. (23) 



Proof. Both identities are obtained thanks to the same kind of considerations. Further- 
more, relation (1221) is obtained in [H Theorem 34] or [23|. We thus focus on identity 

Recall that C2 can be obtained as a limit of 02(77) when rj —>■ 0, where £2(77) is given by: 



C2(r/) : = 



E 



E[Ai(r/) Ai2{ri)] 



We can thus apply identity (1141) with = 1, Si = 0, ti = 1, S2 = 1, ^2 = 2, use expression 
( fT2l) for the kernel /C, and let r/ ^ in order to obtain: 

C2 = \lH ^ |Sl-S2p^^'(s2^-l-|Sl-S2|'^+|s2-l|''') ciS2 rfSi 
'■= C2,l + C2,2 + C2,3 + C2,4, 

with := H{2H — 1). It should be noticed here that, since we are integrating on 
the rectangle [0, 1] x [1,2], the limits as 77 — > can be taken without much care about 
singularities of our kernels [x]^ for negative /3's. Moreover, direct calculations yield 



1 f f fj \ \'? fj 9? 7 f 

C2,l = ^lH J J IS1-S2I ds2dsi = Y J s 



( „2H~l 
1 l-'l 



H 



2 



dsi 



and 

C2,2 ■ 

Finally, we have 



2 rl 



\si — S2p^ ^ ds2 dsi 



C2,3 



and 



1 JQ 



2 



1e 

2 



-7(2 



2H 



1 ^0 



I Si — 52!'^'^ ^ ds2 dsi 



IH 



2{4:H - 1) A ^ 



AH -I 



\si - 1 

2H- 1 
'8{4H-1) 



lAH-l 



2). 



dsi 



(2^^ - 2), 



C2,4 = -IH 

_ H 
~ ~2 



2 i-Y 



1 JO 



S2-\?^\Sx- 82^""-'' dS2 dSi 



2H-2 



S2-ir^(|2-S2 



\2H-1 



-\l-S2 



2i^-l^ 



1 H 



x-l\^^\2-x 



2H~1 



dx. 
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By the substitution y = x — 1 we obtain 

1 Jo 
and moreover, by setting y = —x + 1 we have 

•1 

|2H|r, ™|2_ff-l T„ / „.2Hn , „.l2H-l 







where these two integral expressions appear respectively in the expressions for C2,i and 
C2,4. Hence, putting together our elementary calculations for C2,i, . . . , C2,4, expression (i23l ) 
follows easily. 

□ 

3.2. Proof of Theorem I I.IL Recall that we have reduced our L^-estimates to the evalu- 
ation of where Jj = Ai^i+i. We are now ready to compute those terms, separating 
three different cases: 



(1) Diagonal terms. By stationarity of the increments and thanks to Lemma (3. 11 we 
have 

E[|/.r] = E[\Aoi\'] =ci. 
So (122|) in Lemma O and (EH) give 

g E [\Jlf] = ■ f (mH, 2H) + . n-'^^\ (24) 

(2) Secondary diagonal terms. Using again the stationarity of the increments and 
Lemma [3.11 we obtain 

E = E [Aoi Au] ■ 

Hence by (l23ll in Lemma [3?T] it follows 

n-l 

J2 E [Jfjr^,] = 2T'^ .c,-{n- l)n-'^. (25) 

«J=0,|i-j|=l 



(3) OfT-diagonal terms. Let us consider now the off-diagonal terms, which will induce 
most of the differences in the L^-limit according to the value of the Hurst parameter H. 
Observe first that, as in the proof of Lemma [2?6l for |i — j| > 1 it holds: 

E[Ai,i+i{T])Ajj+i{T])] = / / / K'{r]]Si,S2)K'{ri]Ui,U2)duidu2dsids2. 

J i J j J i J j 

Since we are now away from the diagonal, one can take safely the limit 77 ^ in the 
expression above, which gives: 

E[A^,^+lAj,J+l] (26) 

= H^i2H-lf / / / \ui-U2\'^^-^\sr- S2\'^^-^dUidU2dSidS2. 

J i J j J i J j 

Now we have to distinguish between four cases: 
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\si- S2\^^~'^ dsids2] -n-*^ 

'3 



\AH-A 



(i) The case H < 1/2. Equations ^ and (ES) yield directly 
|E [Jfj;] I < T^^ \^H{2H - 1) jf 

Since |« — j| > 1, the mean value theorem gives 

|E [Jfj;] I <Cn-^^|^-j-ir 
where C is a constant depending only on H and T. Note that for < 1/2 we have 

5^ K-j-ir-^<oo 

l*-il>i 

and so it follows 

J2 lEJfj;! = o(n-^^+^). (27) 

l«-il>i 

('mJ T/je case 1/2 < H < 3/4. Applying again relation fl2Tl) and the mean value theorem 
to the integral on the right hand side of relation (l26l) , we obtain 

|j-j|>l |i-il>i 

where ^ !)• Note now that 

|i~j|>l \i-j\>l l«-il>l 

Moreover, it is readily checked, thanks to a Taylor expansion together with the fact that 

E|,_,|>iK-J-l|'^-'<oo, that 

E K-j±ir-'= E \^-jr-'+o{i). 

\i-j\>i K-il>i 

Hence, we have 

E Efjfj;] =:^i/2(2^_^)2^-4H J2 |z-jf^-^ + 0(n-^^). 
|i-j|>l li-il>i 
Now, observe that 

71" "1 I' "2 2 

E \^ - j^-' = 2 E E 1^ - ^1'''"' = 2 E E^""""' 

ii-i|>l i=2 j=0 1=2 j=2 

n—1 n—l n—l 

■4H-4 



2^^^.4^-4^2E(n-l-j)/ 



n—l n—l n—l 

2nE/^"'-2E^''''^'-2E^"''"' 

j=2 j=2 j=2 

n—l n—l n 

2n E j'"""' - 2 E/"""' - 2 E^''"""' + 0(1)- 
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Let US treat those 3 terms separately: since H > 1/2, we have AH — 3 > —1, and thus, 
by Riemann sums convergence, we get 



n-l 



lim > 

n. — ^riCi f ^ 



n— ^OO 



4H-3 ^1 







4:H-2 

It is thus easily seen that n~^^ X]J=i J^^ '^ — Oin'"^). Moreover, since 

n-l 

where C, stands for the usual Riemann zeta function, we have 

n n 

2n-4^ ^nj^^-^ - 2n-^^ J]/^-^ = 2C(4i7 - 4) ■ n'^^^^ + o{n- 
j=i j=i 
So altogether we obtain 

E [j/^j;] = ^C(4 - 4H)H\2H - 1)^ ■ n'^^+i + o(n-^^+^). (28) 

K-il>i 



-4i^+l^ 



-1/2 



(^Mzj T/ie case = 3/4. Proceeding as in the previous case we obtain 

\i~j\>l K-j|>l 

= i^H\2H-lfn-' Y \^-J\-' + 0{n-'). 

|i-j|>l 

Moreover, following again the computation of our Case (ii) above, we obtain 

rt— 1 n— 1 

J2 \^-J\-' = 2J2i^-l-J)r' + Oil) = 2J2i^-l)r'-2{n-l) + 0{l). (29) 

|i-il>i i=i i=i 

Clearly, 2n~^{n — 1) = 0(n~^). Moreover, since J2]=i = c + log(n) + o(l), where c 
stands for the Euler-Mascheroni constant, we get 

n-l 

2n-3E(n- l)j-^ = 2n-Mog(n) +0(n-2). 

Hence, plugging these two relations into equation (l29l ). it follows 

E lEJf j;| = -^H\2H - If log(n)n-2 + 0{n~^). (30) 

l«-il>i 



T/te case H > 3/4. Along the same lines as in the previous cases, we end up with: 



rp4H 



|i-j|>l |i-il>i 

l*-il>i 
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Since AH — 4 > —1, we now obtain: 



E 

N-il>i 



I- J 



n 



n 



1 /•! 



JG 



\x — y\^^ ^ dx dy 



{AH -3){AH -2)' 



which yields, for any H > 3/4: 

K-il>i 



T^" H\2H-1) 
~A AH-'i 



n ^ + o(n ^) 



(31) 



Theorem O now follows easily from combining (l24ll. (l25l). (l27ll. (l28l). (!30l) and (|3T]). 



3.3. Proof of Theorem II. 2L Recall that the Milstein-type scheme is given by 

^ n— 1 

J o / ^ V iT jn {i+l)T/nJ \ {i+l)T/n iT/n) 



i=0 



As in the proof of Lemma 13.11 using the scaling property, the self-similarity of fBm, 
Lemma 2.6 and moreover letting rj and applying dominated convergence (note that 
we assume here H > 1/2), the mean square error of the Milstein-type scheme is given by 



71—1 n— 1 



EE/ / Oij{si,S2)\si - S2\'^"~'^ ds2dsi 

i=Q j=0 

with -iH = H{2H - 1) and 

^..(^1, S2) = iE(25(:) - 5f ) - fiS)(25W - 5f ) - fij^,) 

for si e + 1], S2 e + 1], i, j = 0, . . . ,n - 1. 

We first show that the contribution of the off-diagonal terms to the error is asymptot- 
ically negligible, i.e.. 



n 



-4H 



2n 



|t-j|>log(n) 



/ 0,j{suS2)\si-S2\^''' 
J dij{si,S2)\Si - S2 



ds2 ds\ 



(32) 



2H-2 



dS2 dS] 



In [17] (see Appendix A) it is shown that 

"i+i i-j+i 



oin 



-AH+l\ 



«-i>log(n.) 



6i,j{si, S2) ds2ds-i 



iH-2 2H-2 



< C ■ (log(n))^^-' ■ n 



(33) 



To use this estimate, define now 

Rlj = {si, S2 e[i,i + l]x + 1] : 9ij{su S2) > 0}, 
{si,S2 e[i,i + l]x [j,j + 1] : eij{si,S2) < 0}. 



R 
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An application of the mean value theorem gives 

+ j ^ 0,j{si,S2)ds2dsi + \i-j-l\^^~^J j ^ eij{si,S2)ds2dsi 

M+l l-i + l 

< / 6i,j{Sl,S2)\Si - S2\'^^~'^ ds2dsi 

Ji Jj 

< Iz-j-lp^-^y g^^^^si,S2)ds2dsi + \i- j + J j ^ Oijisi, S2) ds2ds 

Note that 

Oijisi, S2) ds2 dsi 



J j ^ dij{si,S2) ds2dsi 



< 2, 



so it follows 



n 



-4H 



J 6i,j{Sl,S2)\si - S2\'^"~'^ ds2dsi 



with 



\pi,j\ < C • |i - j - 1 



\2H-3 



Using (l33l) we thus have 



n 



-AH 



/ di,jisi,S2)\Si - 82]'^"''^ ds2dsi 

< C ■ (log(n)e^-^n-2^-^) + ^'"""iP^^^- 



\i-j\>log(n) 



Since 



K-jf''"'<^ ^'''"' = C'(nlog(n) 

l*-il>log(") i>log(n) 

we have obtained 



n 



-AH 



/ Oij{Si,S2)\Si - S2\^^~'^ ds2dsi 

that is (1321) . 

(ii) Now consider the "close to diagonal" terms. Since 

X (l[i,^2](M2) - I[s2,i+l]{u2))\ui - Us\^^~'^ dU2du 

we have for |« — j| > 1 that 

%M^S2)\<C-\l-J-l\^''-\ 



o[n 
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Thus it follows 



l<|i-i|<log(n) 



If if < 3/4 then 



Moreover, if if = 3/4 then 



l<l*-il<log(n) 



4H-4 



<Clog(n)^ 



■AH-4 



i=l 



i=l 



< OO. 



i=l 



where c is again the Euler-Mascheroni constant. Finally, if if > 3/4 we have 

n ^1 



i=l 



1 



4iJ-3 



and so 



Hence we obtain 



■4H-4 



0(n 



4//-3N 



n 



-4H 



/i+l H + 1 
I 6'ij(si,S2)|Si - S2p^"^(iS2(iSi 



l<|i-j|<log(n) 

('mj Combining step (i) and (ii) we have 



o[n 



n 



-AH 



J Oi,j{Sl,S2)\Si- S2\'^^~'^dS2dSi 



o(n 



-4J^+l^ 



l<|'i-i|<n'" 

Therefore, it follows that the leading error term of the Milstein-type scheme is given by 

IhT''' / /. 0,,{s,,S2)\s,-S2\^''-^dS2ds, 

0<|i-j|<l 

n-1 n+l 



"1 1 

ii-ji=i ^ 
Using again the scaling and self-similarity property of fBm we obtain 

t-i+l pj+l 

lim n^^'i y / / 9ij{si,S2)\si-S2\^^~^ds2dsi = ai{H) 
o<|i-il<i * ^ 



DISCRETIZING THE FRACTIONAL LEVY AREA 



17 



with 

which is the assertion of Theorem ll.2[ 



4. Asymptotic error distribution of the Euler scheme: if < 3/4 

Let us first explain the strategy we have adopted in order to obtain our central limit 
theorem for the difference Xt — XJ^ of the Euler scheme and its approximation in the case 
H < 3/4. First, recall that the random variable Xt — X^ can be expressed as 

" l'{i+l)T/n 

Xr-X- = Y,Jr, with Jr ^ / - Blfj dBf) . 

~l JiT/n 

With this expression in hand, it can be seen in particular that Xt — XJj. is still an element 
of the second chaos of our underlying fBm B. 

Let us then recall the following limit theorem for random variables in a fixed finite 
Gaussian chaos, which can be found in [22], Theorem 1]: 

Proposition 4.1. Fixp > 1. Let {Z^, n > 1} be a sequence of centered random variables 
belonging the p^^ chaos of a Gaussian process, and assume that 

hm E[Z2] = L (34) 

n— >oo 

Then converges in distribution to a centered Gaussian random variable if and only if 
the following condition is met: 

lim E[Z^] = 3. (35) 

n— >oo 

This is the criterion we shall adopt in order to get our central limit theorem. The second 
order condition (1341) is simply a normalization step, so that the essential point is to analyze 
the fourth order moments ofXT—X^ in order to prove condition (i35ll . It should be stressed 
at this point that [22l Theorem 1] contains in fact a series of equivalent statements for 
condition (i35D . based either on assumptions on the Malliavin derivatives of the random 
variables Z„, or on purely deterministic criterions concerning the kernels defining the 
multiple integrals under consideration. These alternative criterions yield arguably some 
shorter computations, but we preferred to stick to the fourth order moment for two main 
reasons: (i) The computations we perform in this context are more intuitive, and in a 
sense, easier to follow, (ii) As we shall explain below, the fourth order computations lead 
to some visual representations in terms of graphs, and we will able to show easily that the 
CLT is equivalent to have the sum of the connected diagrams tending to 0. As we shall 
see, this latter criterion is really analogous to [22l Theorem 1, Condition (ii)]. 

In the remainder of this section, we check condition (l35]l for Xt—XJ^, rescaled according 
to Theorem 11.11 in order to get a central limit theorem for our approximation. We shall 
first explain the basics of our diagrammatical method of computation and show how to 
reduce our problem to the analysis of connected diagrams. Then we split our study into 
regular and singular terms. 
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4.1. Reduction of the problem. Owing to Theorem 11.11 it is enough for our purposes 
to show that hm^^oo E[Z^] = 3, where 



n 



(36) 



i=l 



and where the index i varies in {1,2} according to the value of H. Furthermore, the 
self-similarity of fBm implies that 



, j=i 



where Jj = Ai^i+i is the Levy area between i and z + 1. Now, the most naive idea one 
can have in mind is to write Z„ as lim„_,oo ^n l??), where Z„ is obtained by considering 
regularized areas based on -8(77), and then expand E[Z^(?7)] as 



{aAH)n) 



t{H)n) 
4 

n 



ii,...,i4=l b=l 
ij+l r-Xj 

dxj / di/j 



(37) 



E 



.i=i 



E 



where we have used formula ( 1151 ) with = 2 in order to get the last equality. 

We apply now Wick's formula f |T4D in order to get an expression for the expected 
values above, and this is where our diagrammatical representation can be useful. Indeed, 
E [11^=1 iv)] E [11^=1 By) iv)] is the sum of 9 different terms, connecting the Xj's two 
by two according to formula (fTl]) . and also the ?/j's two by two. Each term may be 
represented by a four-point diagram in the following way. Draw a simple line, resp. a 
dashed line between i and j if Xi and Xj, resp. yi and yj are connected. This procedure 
yields 9 different graphs, whose typical examples are given at Figure [U Moreover, the 





34f 



Figure 1. Two examples of diagrams. 



reader can then check easily that diagrams fall into two types: connected ones (6) and 
disconnected ones (3). Furthermore, up to permutations of the indices, there is only one 
disconnected diagram, namely the first diagram of Figure [U One checks immediately 
that the corresponding integral is E[Iij^{ri)Ii^{ri)]'E[Ii^{ri)Ii^{ri)] . Write also the total 
contribution of the 6 connected diagrams as 'E[Ii^{'r])Ii2{ri)Ii.^{7])Ii^{r])](^c)- Thanks to our 
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graphical representation, it is then straightforward to prove the following: for arbitrary 
1, . . . , n, we have 



constants Q, z 



E 



j=i 



3E' 



E 



(38) 



Hence our condition (1351 ) is satisfied for Z„ defined by ( 1361) if and only if the right-hand 
side of equation fl38l) goes to zero for q = n~^/^ (q is in fact independent of i). It should be 
stressed at that point that the latter condition (which is what we call connected diagrams 
go to 0) is an analog of criterion (ii) in [22], Theorem 1], but is obtained here without 
Malliavin calculus tools. This terminology is inspired by the Feynman diagram analysis 
in the context of quantum field theory, see e.g. [T3] . 

Let us set now Zn{ri) = J27=i'^iiv)- With the above considerations in mind, we are 
reduced to show that 



1 X. 

lim lim — ;tE 



0. 



(39) 



This relation will be first proved for H G (1/2,3/4). In that case one may consider 
directly the situation where rj = 0, that is the infinitesimal covariance kernel (x, y) i— > 
H(2H — l)|x — yp^^^, since it is locally integrable. The proof requires only a few lines 
of computations. Each diagram in E[Z^(?7)](c) splits into regular terms - which are also 
well-defined for H < ^ - and singular terms - which diverge when H < ^. As we shall 
see, the bounds given for the non-singular terms also hold true for if < |. Then we shall 
see how to bound the singular terms for arbitrary H by replacing the ill-defined kernel 
H{2H — l)\x — y\'^^~^ with its regularization K'{r]; x, y). This step is of course only needed 
in the case if < |, but computations are equally valid in the whole range H G (1/4, 3/4). 
In other words, the barrier ii = | is largely artificial (the proofs of the two cases are 
actually mixed, and one could also have written a general proof, at the price of some 
more technical calculations). 

Before we enter into these computational details, let us reduce our problem a little 
bit more: recall again that we wish to prove relation (l39|) for Znirj) = J2i=i^iiv)- As 
explained above, we evaluate E[Z^(?7)](c) with 6 different connected diagrams. Let us 
focus on the term, which will be called T, corresponding to the diagram given at Figure [2] 
(the other ones can be treated in a similar manner) . Now, starting from expression ( l37l) , 



2 



Figure 2. Typical connected diagram. 

taking into account the fact that we are considering the particular diagram given at Figure 
[2]and integrating over the internal variables y, we end up with T = -^^^ -^(n,...,^), 
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where (recalling that the kernel K is defined by equation (fTTI) ) 

K'{r]; Xi,X2)K'{r]; X3, X4)Ki^,i3{r]; Xi, x^)Ki^^i^{7]] X2, X4). (40) 

The latter expression yields naturally a notion of regular terms and singular terms: 
split the set of indices {1, . . . , n}^ into Ai U A2, where 

A2 = {{ii, . . . ,^4) I 1 < • • • ,U < n, lii-isl, 1^2-^41 < 1}, Ai = {1, . . .,n}\A2. (41) 

Regular terms, resp. singular terms are those for which \ii — i2\, {is — i^l > 2, resp. 
\ii — ^2! < 1 or 1^3 — i^l < 1. Split accordingly the sets of indices Aj, j = 1,2 into 

reg U ^j,sing, ^nd dcUOtC 

'T'j,reg = hh,-M) and 7;-sing = ^ hh,-M)- (42) 



sing 



(h,...,u)(iAY^ (ii,...,M)eA;' 

It remains to prove that 7J-,reg = o(n^) and 7^-,sing = o(n^), for j = 1,2. These two steps 
will be performed respectively at Section 14.21 and | 



4.2. Regular terms and case H > 1/2. This section is devoted to the study of TJ^eg, 
and also of 7J,sing for > 1/2. In both cases, one is allowed to take limits as ?7 — > 
without much care, by a standard application of the dominated convergence theorem. We 
skip this elementary step, and consider directly our expressions for 77 = 0. 

Let us start by 7^,rcg, which is given by 

^l,rog = ^ / dXi--- dX4K'{xi,X2)K'{x3,X4) 

|ii-i3l,Ki-i2|,|i3-i4l>2-^^i 

Ki^^i^{xi,Xs)Ki^^i^{x2,X4). (43) 

We shall bound this integral by different methods in the cases H e (1/2, 3/4) and H < \: 
Assume first G (1/2,3/4). Whenever |s — — j| < 1, recall from Lemma [2771 that 
Kij{s,t) < |t — sp^~^ if \i — j\ > 2, and s G [i,i + l],t G [j,j + 1]. In particular, the 
quantity \Ki-^^i^{xi,X3)\ in equation (l43l) is bounded by \xi — Xsp^"^. We also obviously 
have \K'{xi,X2)\ < \x2 — Xi\'^^~'^ and \K'{xs,X4)\ < \x4 — Xsp^"^. As a consequence, 

rh+l /-M+l 

|7L,reg| < 2C ^ / dXi--- dX4\x2- Xi\'^^-'^\X4~ Xsl^^-'^ 

|n-j3|,|«l-«2l,|i3-«4|>2 



X Ixi-Xsl^^ ^\Ki„i^{x2,X4) 



Let us undo now the initial scaling by setting tj = xj/n. One gets 

i^i^egi < dt,... f dt4 \t2 - t,r-'\t4 - t,r~' 

Jo Jo 

h - ti\^"~^Kint2\,lnU]int2,nt4). (44) 

Applying Lemma [279l to the above expression ( !44l) and integrating successively with respect 
to ti and ts yields 

|ri,,eg| <n^+'^2i/-2) /" ^^f dt4{l + \t2-t4\^''-')KintMnui{nt2,nt4). (45) 
Jo Jo 



DISCRETIZING THE FRACTIONAL LEVY AREA 



21 



Recall now that \K^nt2\,lnt4,\{n't2, nt4^)\ < min(l, {n\t2 — t^lY^ ^). Hence, one can bound 
this kernel by 1 on [0, 1/n] and by (nt)^^~^ on [1/n, 1], yielding 

/ dUKint,\,[nu\ (nt2, nU) < dt + n^^"' / t^"""^ dt < + n^^'^, (46) 

^0 Jo Jl/n 

and also 

/ dU\t2-U\'''-^K^nt,l,lnui{nt2,nU) < / t'^'-^dt + n'''-'' t^^'-^dt 

Jo Jo Jl/n 

Hence one has found: iTi rogl ^ n + n^^'"^ + n^^^^ < n + n^^~^. In particular, if < 3/4, 
then |71,reg| = o(n^). 

(a) Assume now H < ^. In this case, the integrals we have been manipulating above are 
divergent, so that we will use series arguments instead. Let us observe then that, under 
the same conditions as in the case H e (1/2, 3/4), the bound \Kij^^i.^{xi, Xs)! < \ii —iz 
holds true. We also bound the factor \Ki^^i^[x2,Xi)\ by a constant in order to get 



\2H-2 



1^1. 



regl 




il,h:\h-i3\>2 
«1,«3:|«1-«3|>2 

We now leave to the reader the task of checking, with the same kind of computations, 
that iTi^smgl = 0{n) {provided H > \). 

Turn now to the complementary set of indices, A2: by simply bounding the kernels 
Kij{x,y) by constants in (i43l l. one gets 

rii+l rii+l 

\%,mg\ ~ / dXi--- dX4\K'{xi,X2)\\K'{x3,X4)\ 

|n-«3|,|'i2-«4|<l;|«l-«2|,|«3-«4|>2 *i 



*1,«2:|«1— i2|>2 

Hence |72,reg| = 0{n'^^~'^) = oiji^) when H < 3/4, which is enough for our purposes. 
Finally, provided H > ^, some similar elementary considerations prove that 

\T2,sins\<n(^J^ dxi dx2\K'{xuX2)\^ =0{n), (48) 

where we have used the fact that \ij —ik\= 0(1) for j, /c = 1, . . . , 4 if (ii, . . . , ^4) G A2,smg- 

4.3. Singular terms in the case H < |. Let us reconsider the terms Ti^sing and 72,sing 
in (l42l) . taking now into account the fact that we deal with the regularized kernels 
Xi, X2 ), K'{r];x3,X4) instead oi K'{xi,X2), K'{x^,x^). 

In order to treat all the terms appearing in our sums in a systematic way, let us introduce 
a little of vocabulary: consider any multi-index (ii, . . . ,ip), p > 2 (in our case p = 4). 
We shall say that {ij^, • • • , "^jfc}, ji 7^ • • • 7^ jfc is a maximal contiguity subset of (zi, . . . , ip) 



22 



A. NEUENKIRCH, S. TINDEL AND J. UNTERBERGER 



if (up to a reordering) ij^ — ij^ = . . . = ij^^ — = 1 and ii > ij^ + 2 or < ij-^ — 2 if 
/ 7^ ji, . . . , jfc- Maximal contiguity subsets define a partition of the set {ii, . . . ,ip}. Then 
we shall write (ii, . . . ,ip) G Jmi,...,mq if the lengths of the maximal contiguity subsets of 
{ii, . . . ,ip) are nii > . . . > rriq. 

This terminology will help us classify the terms in T^^sing U 7^,sing- Forgetting about 
the 0(n) multi-indices (ii,...,i4) in J4 appearing in 72,sing (according to the fact that 
Vai^Astiv)) is uniformly bounded on [0,T], proved in [25], this term contributes only 
0{n) to the sum), the other singular terms are all in Ti smg and may be: 

- either of type ^2,1,1, with maximal contiguity subsets {{21,22}, {^s}, {^4}} or equiva- 
lently {{23,^4}, {«i}, {^2}}; 

- or of type J2,2, with maximal contiguity subsets {{21,22)5 {^35^4}}; 

- or of type J^^i, with maximal contiguity subsets {{21,22,23}, {24}} or equivalent pos- 
sibilities. 

Let us observe that, in our iterated multiple integrals, the most serious problems of 
singularity appear when the external variables x (represented by solid lines in our graphs) 
are contiguous. Indeed, the internal variables y are integrated, smoothing the kernels 
K' into Ka^i,. However, one still has to cope with the highly singular kernel K' for the 
external variables. For instance, for the graph given at Figure [2] (which is the one we are 
analyzing), this kind of problem appear for the terms of type ^2,1,1 (when the maximal 
contiguity subset is {{21, 22}, {^3}, {^4}}) or J2,2- But a simple Fubini type argument allows 
us to get rid of these singularities. Indeed, when r] > 0, the integral 

[[ / dxj K'{r];xi,X2)K'{7];x3,X4) ■ [[ / dyj Js:'(r/; 2/1, 2/3)K'(r/; 2/2, 1/4), 
j=i ''^j j=i •'^j 

corresponding to the diagram of Figure [2|, is also equal to 

W I dyj K'{T];yuy3)K'{r];y2,yA) ■ JJ / dxj K'{r];xi,X2)K'{r];x3,X4), 

j=l -^ij j=l -^Vj 

corresponding (up to time-reversal) to the reversed diagram obtained by exchanging full 
lines with dashed lines. The important point is that this full-line dashed-line symmetry 
maps the above singular diagrams of type ^2,1,1 or J2,2 into regular diagrams, for which 
the external variables are separated. This situation can thus be handled along the same 
lines as in Section [42l and there only remains to estimate singular diagrams of type J3 1. 

For this latter class of diagram, assume for instance (without loss of generality) that 
{^1) ^2, is} is a maximal contiguity subset of (zi, . . . , ^4). Then, owing to relation (fT2l) . the 
corresponding integral writes E = E{ii, . . . ,14), with 

/■«3+l fii+i /■«2+l 

E = Ch dX3 / dXi / dx2 / dx4 [X3 — X4]'^^~'^[xi — X2]'^^~'^ 

J 13 J ii J 12 J ii 

{[X3 - Xi]f + [Z3 - - [X3 - - [Xl - l3]f) Ki,,iMx2,X4), (49) 

which is the sum of 4 terms, denoted in the sequel by Ei, . . . , E4. The most complicated 
one is a priori £'1, obtained by choosing the contribution of [X3 — xi]^^ to the integral. 
Let us first estimate this term. 

Apply Lemma 12^81 with f{x4]u) = [u — X4]^^~^, z = Xi {x^ is simply an additional 
parameter here, and / fulfills the analytic assumptions of Lemma 12.81 because 13 and 14 
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are not contiguous) and /5 = 2iJ, 7 = 0: letting 

^2H\^ ^ 12H-2 



(j)i{xi; xi) := / dxslxi - Xs].^ [0:3 - ^4]^ 

we obtain that 0i is analytic in xi on a cut neighborhood fi'^^j of + 1] excluding 
possibly is and is + 1 (depending on whether i^, i^ + l & {ii, ii + 1} or not), and one can 
decompose 01 into 

0i(x4; xi) = [xi - i3]J^+^Fi(x4; Xi) + G'i(a;4; Xi) (50) 

on a neighborhood of is (and similarly around is + l), with Fi possibly zero. The functions 
4>i\n'^^^, Fi and Gi are analytic and bounded by a constant times jis — Z4p^~^. 

Apply once again Lemma [23] with /(X4; u) = 0i(x4; u), z = X2 and (3 = 2H — 2, 7 = 
or (possibly) 2H + 1: letting 

pii+l 

</)2(x4;X2) = / dxi [X2 - Xi]J^"^(/)l(x4;Xi), (51) 

02 is analytic in X2 on a cut neighborhood of [^2,^2 + 1] excluding possibly ii and 
ii + 1, and 

02(2:4; X2) = [x2 - ii]l"~^F2{x4,; X2) + [x2 - ii]t," Fs{x4] X2) + G2{x4; X2) (52) 

on a neighborhood of ii (and similarly around ii + 1), with the same bounds as before for 
02b^'„,, F2, Fs and G2. 

Finally, since 02 is integrable with respect to X2 on [22,^2 + 1] and Ki^^i^^rj; X2, X4) is 
bounded by Cjis — i4p^~^ by Lemma [2?7l one gets 

1^1 < G' dxi Ks - lA^"'^ = C'\iz - 24|'^"'. (53) 

J i4 

There remain 3 'boundary' terms E2, E3, E^ which are easier to cope with. Consider 
for instance E^ defined as 



Es = / dxi i dx2Ki2,u{ri; X2,X4) 



x/ dxi[x2-xi]l^~^ dxs [xs-zi]5^[xs-X4]5^~^ 



Applying again Lemma [2.81 we get 

E^ = 

G / dx4Gi{x4] ii) / dx2Ki^,u{'r]; X2, X4) ([xs - H - - [^2 - ii]J^~^) , 
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are 



where Gi is as in eq. ( !50l) . Since X2 ^ [x2 — ii — 1]^ ^ and X2 1— > [x2 — ii]^ 
integrable and Gi, resp. -^'12,44 is bounded by a constant times [is —24^^"^, one easily 
gets an upper bound as the same form as before, namely, |i?3| < Cjis — i^l'^^'*. 

We have thus proved that E{ii, . . . ,14) defined by ( !49l ) satisfies E{ii, . . . ,14) < Cjis — 
^^|4H-4_ Finally, since ^|j3_j4|>2 — 14]^^''^ = 0{n) (as in eq. (!47l)). we obtain 
En,...,^4GJ3,l^(^l'•••'^4) = 0(n). 

Let us summarize now the results we have obtained so far: we have shown, respectively 
at Section [42] and [4731 that the terms TJ.rcg and TJ.smg defined by equation ( [42] ) are o{'n?). 
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Going back to the definition of T (see equation ([40l) ). this also shows that this quantity is 
of order o(r;,^). Recall now that E[Z^(?7)](c) can be decomposed into 6 terms, corresponding 
to our connected diagrams, each of the same kind as the particular example T we have 
chosen. We have thus proved that E[Z^(?7)](c) = o{n^) uniformly in 77, which yields relation 
(l39l) . This finishes the proof of Theorem II .31 for H < 3/4. 

5. Asymptotic error distribution of the Euler scheme: H > 3/4 

In this case, we derive the limit distribution in a different way, and first analyze the dif- 
ference between the Euler and the Milstein scheme. An exact expression for this difference 
is given by 

^ n— 1 

^, -5(f), (54) 

2 / {t+l)/n t/ny^ (i+l)/n «/■«'" V" / 

i=0 

and we will see that, thanks to a simple geometric trick (borrowed from [20]), the latter 
quantity has the same law as 

^ / ^ y («+l)/n «/nl I (i+l)/n i/nl 

1=1 



This allows to apply easily Theorem 2 in [24], yielding the Lemma below, in which the 
following distribution appears: 

Definition 5.1 (Rosenblatt random variable). A standard Rosenblatt random variable 
with parameter Hq = 2H — 1 is given by 

I I ' I [u,s) — — [u, r) du dWrdWs 



4H {2H -if^ Jo Jo Vimaxir.s} du ' du 

where W is a standard Brownian motion, 



Knit, s) = CHs'/'~" j\u - sf~'V~'/' du l[o,)(s) 



and 

Ch 



H(2H -1) ^ 



f3{2-2H,H -1/2) 



Lemma 5.2. The following limits in law hold true: 
(i) Let H = 3/4. Then we have 



^/ci{H) \ogn 

where ci{H) = 9/16 and Z is a standard normal random variable, 
(a) Let H E (3/4,1). Then 

Wiw. - 4';:)(«<?i../» - ^i;^) ^ ;^(«' - «^)- 

where C2{H) = 2H^ {2H — 1) / {4H — 3) and Ri and R2 are two independent standard 
Rosenblatt variables of index 2H — 1. 
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Proof, (i) Let /5 be a fractional Brownian motion with Hurst index H and define 

- ~ I 1 I - -1 + n 2^ |/^(*+l)/n - A/nl • 

i=l ^ ^ i=l 

If = 3/4 it follows from [24] that 



Ci[H) log(n) 

where Z is a standard normal random variable. Moreover, for -ff G (3/4, 1) it is shown in 
1241 that 



^K. ^ fl. (56) 

where i? is a standard Rosenblatt random variable with index 2H — 1. 

Now let (3 be another fractional Brownian motion with the same Hurst index as j3, but 
independent of /9 and define 

n 
i=l 

The continuous mapping theorem and (l55l) implies that 



-— -K -^Z,-Z, (57) 
ci(if) log(n) 

for if = 3/4, where Zi and Z2 are two independent standard normal random variables. 



From f l56l) we obtain that 

^K^(^i-^2), (58) 

where -Ri and are two independent standard Rosenblatt random variables with index 
2/7-1. 

(a) Now, set 5(1) = (/? + ^)/y2 and E^^) = (/? - ^)/v^. Then B^^) and 5(2) are two 
independent fractional Brownian motions with the same Hurst parameter. Moreover, we 
have 

n-1 



Z^'^-"(fc+l)/n ^k/n>\^(k+l)/n ^k/n>~r)^n- 



k=0 



Thus, we have for = 3/4 that 



n-1 



c,{H)\ogn hi^'^'^'-" ^VnK^(.+i)/n ^./n^ " ^ ^^(i/) log(n) 



and the first claim follows from (jSTIl and the fact that Zi — Z2 has the same distribution 
as \/2Zi. 

Moreover, since 

o 2—2H 

\^rR(l) _ r(1)vr(2) _ r(2)n ^ j;?^ T/' 

Z^l-°(i+l)/n ^i/n)y^{i+l)/n ^i/n) ~ / 7TT\^n 



v/^^t^ (^+l)/n Vn/V (.+l)/n ^/n^ ^y^^ n 
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the second claim follows from (l58|) . 

□ 

Since the Milstein scheme has a better convergence rate than the Euler scheme for 
H >3/4:, the error of the latter scheme is dominated by ( l54l ). Thus, the asymptotic error 
distribution of the Euler scheme can be determined by the above Lemma, which will be 
carried out in the following two subsections. 

5.1. Error distribution of the Euler scheme for H = 3/4. By scaling we can assume 
without loss of generality that T = 1. Recall that here we have 

E|Xi - = ^ ■ log(n)n-2 + o{\og{n)n-') . 

for the error of the Euler scheme. Using the Milstein-type approximation X" we can write 
Xi - Xf = Xi - Xf + Xf - Xf 

^ n— 1 

~ n Z_^\^{i+l)/n ^i/n)y^{i+l)/n ^i/nl^Pn, 



i=0 



where p„ = X" — X". Hence, setting := n[-^ log(n)] we obtain 



n-l 



1=0 

Now note that Knpn ^ in L'^{VL) by Theorem 11.21 and 



\^rR(^) r(^)vr(2) R(2)^ ^ 7 



v/ci(i^) log 



1=0 



where ci{H) = 9/16 by Lemma [521 Since [2/ci(/7)]i/2 ^ 1 [128/9] it finally follows 
that 

n{log{n)r'/\X,-Xr)^^-Z, 
where Z is a standard normal random variable. 

5.2. Error distribution of the Euler scheme for H > 3/4. Here we have 

E|Xi - Xi"|2 = ai{H) ■ + o{n~^) 

with 

li/2(2iJ- 1) 



4 AH -3 

Proceeding as above, the limit distribution of the error of the Euler scheme is determined 
by the limit distribution of 



2^fouJ^W)^ '•/"'^^ («+l)/n ilnl 



Since 

n AH -3 V2 



2^^1(17 H^{2H-1) ,/^^- 
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it follows by Lemma [52] that 
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